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1. Introduction 



Conformal invariance serves as a powerful principle in two dimensions. It is tempting 
to generalize what we have learned in two dimensions to higher ones. Some have been done 
along this line 0. Recently, Cardy proposed generalizations of modular invariance to 
higher dimensions f^. He demonstrated modular invariance for a free scalar in various 
dimensions, with geometry S'^xS^. Unfortunately, the form of modular invariance working 
for an anti-periodic free scalar in three dimensions fails in the 3d Ising model at the critical 
point. In this letter, we shall find that this form fails too for a Chern-Simons matter field 
theory in three dimensions. This failure, however, has not yet ruled out the possibility 
that there exist other forms of modular invariance working for this particular conformal 
field theory. A universal expression is even more desirable. 

The quantized Chern-Simons, Abelian and non-Abelian, gauge theories coupling to 
matter fields are believed to be conformally invariant, provided the corresponding classical 
actions involve no any dimensional parameters, including matter masses. It has been shown 
in 1^ that the beta functions of gauge couplings vanish (see as well for the Abelian 
cases) and massless matters remain massless after quantization. Furhtermore, the induced 
non-local term for the gauge field by integrating over matter fields. 



is invariant under the special conformal transformations, x — > x' = i_|_2^x+b2x2 • 

Being conformally invariant makes it possibable to explore the models' other nice 
features such as the modular invariance. For this purpose, let us define the theory with 
Abelian Chern-Simons coupled to massless charged scalar on the geometry S"^ x S^. A 
natural metric on it is 



where dQ^ is the unit round metric on S'^, and L the radius of S^. This metric, viewed 
as a metric on S'^ x R, has a map to the fiat metric in 3d Euclidean space under a Weyl 
rescaling. As the model is conformally invariant, its Hamiltonian on S'^ has a map to the 
dilatation operator on R^. Then, the problem of calculating the partition function of such 
a system on 5"^ x S"^ is equivalent to that of counting the scaling dimensions in the model. 
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ds^ = R^dVt^ + L'^dT^, 
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Our main result in this paper is the foUowing partition function for the model of interest, 
in the infinite k (the statistical parameter) limit 

1 / °° 1 

^^^^^2^iJ^ n (1 _ ^qn+l/2^2n+l(^l _ ^- 1 gn+1/2) 2n+l ' 

where q = exp{—27r ^) is the modular parameter. The contour of the above integral is the 
unit circle on the complex plane. It is an interesting mathematics problem to evaluate this 
contour integral. It should involve certain analog of Rogers-Ramanujan identities. The 
evaluation of this integral is necessary in finding a possible form of modular invar iance of 
the model. 

For the Chern-Simons matter theories with finite k, expression of the partition function 
seems to be much more complicated, because the scaling dimensions of operators may 
not be the canonical ones. It has been seen in Q that the scalar field has indeed an 
anomalous dimension, proportional to -p-, while the Abelian Chern-Simons gauge field 
and so the statistics parameter k need no infinite renormalization. Renormalization of 
some composite gauge invariant operators with a dimension not larger than 3 has been 
investigated in a recent work . It is remarkable that the anomalous dimensions of these 
composite operators have simple relations with that of the primary field. This perhaps can 
be attributed to that Chern-Simons matter theories have vanishing beta functions. If these 
simple relations remain for all gauge invariant operators, a similar result, of k dependent, 
could be obtained. Investigations along this line are in progress. 

The plan of this paper is as follows. In the next section, we shall compute the partition 
function by counting scaling operators. Modular invariance proposed by Cardy is numer- 
ically disproven. In section 3, we calculate the partition function in the infinite k limit, 
upon regularizing a determinant. The result is compared to that obtained by counting of 
scaling dimensions in section 2. The agreement of the results indicates the consistency of 
the approach used. 

2. Partition Function on the Geometry S"^ x 

Just as in a two dimensional conformal field theory, the partition function of a higher 
dimensional conformally invariant field theory on a certain geometry encodes information 
of scaling operators. In two dimensions, one considers usually a torus. The partition 
function on a torus, modular certain conformal anomaly related terms, can be considered 
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as the generating function of counting of scaling operators. This can be generahzed to 
higher dimensions with the geometry S'^ x S^. The natural metric on 5''^ x 5"^ is the one 
introduced in the introduction, provided dQ^ is replaced by the unit round metric dO^ on 
S'^. Such a metric can serve as a metric as well on the cylinder S'^xR. The latter geometry, 
under a Weyl rescaling, is related to the flat geometry in the Euclidean space with the 
origm removed, i?'^+^\{0}. To see this, make a change of coordinates r = (i?/L)logr, the 
metric becomes 

ds^ = ^{dr^ + r^dnl) 

which is conformal to the flat metric on the Euclidean space. The Hamiltonian H on space 
5"^ is then proportional to the dilatation operator D in the d + 1 dimensional flat space, 
up to a conformal anomaly. The partition function on 5"^ x 5"^ is the generating function 
of scaling dimensions x of the model: 

^ (2.1) 

q =exp{—2nd) = exp(— 27r — ), 

R 

where the sum is over all scaling operators. In a recent interesting paper , generalizations 
of modular invariance to higher dimensions were made. The proposal of Cardy works 
for a free scalar field coupled conformally to the scalar curvature of the geometry. The 
particular form of modular invariance, which works for an anti-periodic free scalar in three 
dimensions, does not work for the 3d Ising model at the critical point. Here we shall 
see, to our disappointment, that this form does not work either for the model of Abelian 
Chern-Simons coupled to a massless complex scalar. This, though, has not ruled out yet 
a universal form of modular invariance for 3d conformal field theories. As we do not 
understand a possible underlying principle of modular invariance in higher dimensions, it 
is conceivable that a universal form is far from being simple. 

Note, that the formula ( |2.1| ) works only for a strict conformal field theory. It is 
then a nontrivial check of the conformal invariance to calculate the partition function by 
independent means. We shall calculate the partition function of our model in this section, 
starting from ( |2.1| ). If our conjecture of finiteness is true, then this partition function is 
the one for all coupling k and independent of k. In the next section, we shall calculate the 
partition function in the infinite k limit, starting with a path integral. We shall see the 
results are same. 
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The Euchdean action of our working model is 

S = Jd^x (^iD'(P)+{D,(P) + i^e'^^A.djA^ , (2.2) 

where = 9* + lA^ . The scahng operators in a gauge theory are aU gauge invariant 
operators. However, not all these operators in the model (|2.2D are independent. The 
inter-relationship is specified by the equations of motion 

Ae^^'^F,, = ((!)»+./> ^^^^^ 

These equations of motion will simplify counting of gauge invariant scaling operators. 
The first equation in (|2.3|) means that the strength of the gauge field is not independent. 



So in constructing a gauge invariant operator we need not consider F and its covariant 
derivatives. The second equation in (|2.3|) implies that not all three {Di)^(p (and their 
conjugates) are independent. We can forget, for instance, (il>3)^(^ in our construction of 
scaling operators. 

Now, although Di and Dj with i ^ j do not commute, DiDj is not independent of 
DjDi. This follows from DjDi = DiDj + [Dj,Di\ and the second term on the r.h.s. is 
proportional to Fij which have been given away from the counting of operators, by use of 



the first equation in (|2.3|) . Once all these have been taken into account, a general gauge 



invariant scaling operator takes the form 



/(I) /(I) /{I) , „(i) „(i) „(i) \ / /C"^) 7{™) /(") , 

'p[^ (p)+D^^ D^^ (/> ••• (^i' D2 Di (p)+D^^ D^^ D 
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(2.4) 

While if, if {n{ \ nf) are arbitrary non-negative integers, lf{nf) = 0, 1 is restricted. 



because of the second equation in (|2.3|). Let us assume that the scaling dimension of the 
above operator coincides with its canonical one, as it is when k goes to infinity. The scaling 
dimension of the operator (p. 4]) is 



a: = m-f^(/f)+n!"): 



i.a 



An observation of (|2.4| ) suggests that the partition function of interest can be written 
as a sum 

Z{q) = J25iJ2n - 5^.,)5E.-(V2+n«+n«+n«)+E.-(l/2+4"+4^^+4^\ (2.5) 



where and Si are multiplicities. The introduction of a delta function is due to the fact 
that a gauge invariant operator must contain the same number of (p and (p'^ . This delta 
function renders a direct calculation of the partition function difficult. Let us introduce 
another parameter z, and define 

Z(a z)=^zl^ '^-T. n(l/2+n«+nW+nW)+^ (1/2+Z W +Z« +Z W ) 

^ ^ (2.6) 

=Zi{q,z)Zi{q,z~^). 

where 

Z,{q,z) = 5^^E-gE-(V2+r.«+n«+n«). 

It is easy to see that Z{q) is just the coefficient of the term z^ in Z{q^ z). Now Zi{q, z) is 
easy to calculate, the result is 

1 oo ^ 

Zl{q, ^) — n 1 _ _2Ql/2+ni+n2+n3 ~ 11 (1 - 2Qn+l/2 -) 2n+l ' 
ria n=0 \ ^ ' 

It is interesting to compare this with the partition function of a free scalar, obtained by 
Cardy in §, 

oo ^ 
Zo{q) = n (l_on+l/2)2n+l' 

we see that ZQ^q) = Zi{q, z = 1). 

Finally, the partition function for the model (|2.2|) (at least in the infinite k limit) is 

n=0 ^ ^ ' ^ / 

The contour integral in ( |2.7|) is formally defined. We shall duplicate the formula in the 
next section by working out the path integral in the infinite k limit. We shall see that 
parameter z plays a role of the holonomy of a fiat gauge field along the Euclidean time 
direction. (|2.7| ) is the starting point to understand modular and other properties of the 
the theory. An exact evaluation of it is on the way. 

Now we turn to the numerical investigation of modular invariance in the form proposed 
in . Our conclusion is that this form of modular invariance does not work for the theory. 
As the first step, we check modular invariance associated with a partition function 



oo ^ 

^f"' = n 7131 



n=l ^ ^ ' 



Z{q) is not invariant under 5-^1/5. Write 



logZ(g) 



/ 

Jc 



C(2+e) 



S-'F{s)ds, 



where a factor l/{2ni) is absorbed into the definition of the contour integral. The above 
formula can be obtained by using the inverse Mellin transform of the gamma function, as we 
shall do for the integrand in (|2.7|) in the next section. The contour C(2 + e) is the one with 
Res = 2 + e, and e is a small positive number. The above contour is chosen such that all 
poles of the integrand are located to its left. Explicitly, F{s) = {2Tv)~^r{s)C{s — l)({s + l). 
C{s) is the Riemann zeta function. Now, the quantity 



is modular invariant, except a few additional terms in the form 5", a = ±3/2, ±1/2, will 
appear after the modular transformation 6 1/6. These terms appear due to the poles 
of the integrand in the above contour integral at s = ±3/2, ±1/2. 

Note that, because of appearance of terms 5^^/^, the function $3 defined in an ap- 
pendix in is not modular invariant. $3 = {d^ — 9/4)1^ is defined such that additional 
terms 5^"^/^ are removed, here x = logS. Instead, we define a new function 



This function is invariant under x — >■ —x, thus an even function of a;. If we have the 
following expansion 




(2.8) 



(2.9) 



logZ(g) = 5^aA?\ 



A 



then the function can be expanded as 



A 



To numerically check modular invariance, we expand Z{q) up to g^: 



Z{q) - 1 + g ± 3g2 + + 13g^ + . . . , 



and to the same order 



^ ^ 5 2 10 3 21 4 

logZ(g) =q+-q^ + —q^ + — g^ + . . . . 
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Because the modified Bessel function Ki{z) decays exponentially with the growth of its 
argument, it is sufficient for us to take the first few terms. We used Mathematica to 
calculate function The result is 



= 0.3358 + 5.8 x lO'^x - 1.7267x2 + 6.2 x 10"^a;^ + 4.2673a;^ + 1.6 x 10"^a;^ .... 

It is seen that to a good approximation, is an even function. To see how inclusion of 
higher order terms improves the situation, we expanded Z{q) to and calculated 
We found 

$(a;) = 0.3358 + 1.4 x lO'^x - 1.7267^^ + 3.4 x lO'^x^ + 4.2676x^ + 2.2 x lO'^x^ + .... 

This improves the vanishing of the coefficients of odd terms by a factor 10~^. 

We do not know a model with a partition function investigated above. Cardy consid- 
ered the partition function of an antiperiodic scalar. The antiperiodic condition is crucial 
for modular invariance in that case. Here in our model of Chern-Simons coupled to a 
massless complex scalar, we need not consider such boundary condition, since all gauge 
invariant operators given by ( |2.4| ) is even under this boundary condition. For this model, 
the candidate of a modular invariant function is $ defined above, while not $3 introduced 
by Cardy. The reason is the following. First notice that to fix the additional finite terms 
under modular transformation, we need pick out poles of the integrand in (|2.8|) between 



(7(3/2 + e) and C(— 3/2 — e) 0. It is readily seen that one pole is s = 1/2, provided it is not 
cancelled by a zero of F{s + 1/2). It happens that for a free scalar, this pole is cancelled. 
Then one looks for a possible pole of F(s + l/2) at s = —1/2. Precisely because of the 
anti-periodic condition, this pole does not occur. Now for our toy model discussed above, 
it occurs. If one lift the anti-periodic condition, just as in the model under consideration, 
one normally expects that this pole must be there. Modular invariance requires that the 
integrand (excluding d~^) in (|2.8D is invariant under s — s . So one expects another 
pole at s = 1/2. Other poles at s = ±3/2 will normally occur. We do not expect poles at 
integers. 

Even when there are no poles s = ±1/2, it is safe to use This is because if $3 
defined in is invariant, $ is invariant too. The partition function (|2.7| ), up to order q^, 
is 

Z{q) = l + q + 7q'^ + 26q^ + 82g^ + 233g^ + 657/ .... (2.10) 
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We first use the expansion up to in calculating $(a;). The result is 

= 0.4094 - 0.5971X + 0.2459^^ - 3.0990x^ + 5.7426x^ + 1.9280x^ + . . . . (2.11) 

The coefficients of odd terms are of the same magnitude as the ones of even terms. This 
clearly rules out modular invariance. Next we use ( |2.1C| ) up to order to calculate $(x). 
We find 

$(x) = 0.4094 - 0.5971X + 0.2460x2 - 3.0994x^ + 5.7453x^ + 1.9163x^ .... 

So inclusion of higher order terms merely improves higher terms in $(x), within our ap- 
proximation. The term x, say, is stable when more and more higher terms are included. 

3. Regularizing Determinant 

In this section, we shall compute the partition function of our model by regularizing 
the determinant arising from the path integral in the infinite k limit. To compare the 
result in this section with ( |2.7| ), we shall first recast the logarithm of the integrand in ( |2.7| ) 
into a contour integral. We notice first that 

\ogZ{q,z)= Yl (2n + l)-(2'^ + z-'')Q''("+^/2)_ 

n=0,r=l 

Using the inverse Mellin transform of the gamma function 

Jc 

where the contour C is parallel to the imaginary axis and right to it, and again a factor 
l/(27rz) is absorbed into the definition of the contour integral, we find 

\ogZ{q,z)= [ 6-'{2n)-'{2' -2)T{s)C{s-l){F{z,s + l)+c.c.)ds, (3.1) 

ic(2+e) 

where ({s — 1) is the Riemann zeta function. The generalized zeta function F{z, s + 1) is 

^(^'- + 1) = E;^- (3-2) 

n=l 

The contour in ( p.l| ) has been chosen so that all poles of the integrand are located to its 
left. 
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We now calculate the partition function with the path integral. In the large k limit, 
only flat gauge field configurations contribute to the path integral. In the spacetime with 
topology S'^ X S^, these configurations are parametrized by an angle variable 9 G [0,27r]: 

Note, that we have assumed that the radius of S"^ is 1 so the parameter 6 = L/R = L. 
The circumference of is 271 L. We will see that the complex parameter z introduced 
previously is related to 9 through z = exp{i9). 

Now the partition function in the large k limit reads: 

Z{q)= [ ^det-'/^[-id + iAi9))id + iA{9))+^Rl (3.3) 
J 2n 

where we assumed that the scalar is conformally coupled to the scalar curvature, therefore, 
^ = 1/8 [H. The eigenvalues of A on the unit two sphere are /(/ + 1), / being non- negative 
integer. Formally, the partition function takes the form 



-exp - + V2)log((/ + 1/2)2 ^ ^-2(^ ^ ^/(27r))^ 



(3.4) 



l=0,n= — oo 



In it, the boundary condition along the time direction is not important, since it can always 
be put into the 9 parameter. In other words, one can say that the holonomy of the flat 
gauge field shifts the boundary condition. 

To be well-defined, ( ]3.4D needs to be regularized. We use the ^-function regularization 
to do so for the sum. Define 

°° 1 + 1/2 

The regularized partition function is 

Using formula 

f oo 

-Xt jS — 1 



/•OO 

x-'T{s) = / e-"V 
Jo 



dt. 



we rewrite the zeta function as 



oo 



as,9)r{s)= I e-^F^{5-H)F2{t)dt, (3.5) 
9 







where 

oo 



n= — oo 
oo 



2' 

1=0 



({s,9), defined as above, has obviously a pole at s = 3/2. On the other hand, s = 
is a regular point of (■(■§, d)- Actually, C(0, 9) = 0, since the r.h.s. of ( |3?5| ) is regular while 
r(s) is divergent at s = 0. Denote the r.h.s. of (|3.5|) by g{s). Then 



V{s) T{s) 

We shall see that (/'(O) is regular, so that the first term on the r.h.s. of the above equation 
is zero at s = 0. At s = 0, the second term is 

-g{s) ^ ^ (7(0). 

We thus find 

C'(0,^)=^7(0). 

This is what we need to calculate. Thus, we need to extend the function g{s) in (|3.5|) to 
s = 0. When t — > 0, both Fi{t) are singular. To extract the singular behaviors of these 
functions, one can use integrals to approximate sums: 



Flit) / dxe-"" ' ^ 

Jo Vi 

pQO -I 

F2(t) ^ / dxxe"^'* -. 
Jo t 

To go beyond this approximation and actually calculate Fi{t) for small i, we apply the 
inverse Mellin transform and obtain. 



Fi{t)= / t-'V{s){C{2s,u) + C{2s,l-u))ds 

'C(l/2+e) 



t-^'^V{s/2){as,u)+asA-u))ds, 

C{l+e) 



(3.6) 



F2{t) = / t-'{2''-' - l)r(s)C(2s - l)ds 

Jc{l+e) 

= \ I t-^'\r-^-l)T{s/2)as-l)ds, 
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where u = 9/{2n), u) is the generaUzed Riemann zeta function ||^. The above contours 
are chosen such that aU poles of integrands are located to the left of the contours. 
We then consider t 1/t: 



Fiil/t)=l [ t'/^Tis/2)ias,u) + asA-u))ds 

=Til/2)Vi+l [ t^/^T{s/2){as,u) + asA-u))ds (3.7) 

^ Jc(-l-e) 

=v^t+l [ t-'/^r{-s/2)ia-s,u)+a-sA-u))ds. 

Denote the second term in ( |3.7| ) by It is easy to see that when t gets large, 

decreases exponentially. 
Similarly, we have 

^ ' ' 2 24 960 

+ I f t-^/\-'-y'{2-^-' - l) p/^7/^l,, r(./2 + 1)C(. + 2)ds. 

It is easy to show that the last term behaves like when t gets large. Denote this 
function by f2{t)- 

We are ready to extend g{s). It is 



oo 

g{s)= I t'-^F,{6-h)F2{t)dt 
Jo 

= J t'-^Fi{5-h)F2it)dt+ I t'-^Fi{6-^t)F2{t)dt 

{t-^-'[sv^t + his'm- + -- ^t-' + f2{t)]+t^-'F^{6-h)F2mdt 



--V7tS{- + 



2s -3 12(2s-l) 480(2s+l)' 
, f 1 7 

+ {t-^-'[sV^tMt) + C- + -- ^t-' + Mt))hi5H)]+t^-'F^i5-H^^^ 

(3.9) 

And 

207 

/oo / 1 7 (3.10) 

t-'[F,{6-'t)F2{t) + 6V^tMt) + (2 + ^ - ^t~' + f2{t))h{6't)]dt. 
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According to the properties of functions fi{t) at large t mentioned above, the integral on 
the r.h.s. of ( |3.10| ) is well defined. We shall evaluate this integral term by term. First, we 
calculate 



/oo 
t-^Fi{5-^t)F2it)dt = ay/iiS 

+ / ds5-'{2'-^ - l)(27r)-T(s)C(s - l)[F{z, s + l) + c.c] 

Jc(2+e) 



(3.11) 



where 



/ -^{r-'-ms/2)as-i)ds. 

Jc(2+e) 'S + J- 



'C(2+e) 

We have used Hurwitz formula in evaluating (|3.11|) 

oo 

C(-s, u) = 2(27r)-"-^r(s + 1) XI ^"^"^sin(2™ - ^s). 

n=l 

For the first time we have seen F{z, s + 1) introduced before in ( |3.2| ), z = exp(i9). Other 
terms can be calculated similarly. The final result is 

C'(0, 6)= [ 5-^(2^ - 2)(27r)-^r(s)C(s - l){F{z, s + 1) + c.c). (3.12) 

Jc{2+e) 

This formula is exactly the same as in ( |3.1|) . 
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